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Abstract. The changes in the spectral line profile of the 22 GHz H2O maser are calculated as a function of emerging maser 
flux. We address not only the narrowing and re-broadening of the maser lines, but also the deviations of Gaussian symmetry as 
a result of the various hyperfine components of the 22 GHz maser line. Non-UTE models of the HtO maser transition, including 
hyperfine structure and cross-relaxation, are compared with high spectral resolution Very Long Baseline Interferometry (VLBI) 
observations of the H2O masers in the circumstellar envelopes (CSEs) of a sample of late-type stars. This yields estimates on 
the thermal width in the maser region as well as the emerging maser flux and thus the level of saturation. We also discuss the 
effect of a velocity gradient along the maser path on the line widths and shapes of the line profile and the effect of the geometry 
of the maser region. We find that the typical velocity shift along the maser path is of the order of 1.0 km s"'. The effect of this 
shift on the shape of the maser spectrum is difficult to distinguish from the effect of the hyperfine components. 

Key words, masers - stars: circumstellar matter - stars: individual (U Her, S Per, NML Cyg, VY CMa) - stars: AGB and 
post-AGB - Line: profiles 



1. Introduction 

The spectral line profiles of astrophysical masers contain im- 
portant information on the conditions of the masing region, 
such as the thermal line width of the particle distribution and 
the optical depth. However, the degree of saturation of the 
maser has a strong influence on the maser profile and can- 
I not be determined directly. It has been shown before, that a 
detailed analysis of the maser line widths can provide some 
information on the maser saturation (e.g. Goldreich & Kwan 
1974; Litvak 1970). The line profiles first get narrower with 
increasing radiative flux and then re-broaden when the maser 
becomes radiatively saturated. However, if there is an effec- 
tive redistribution of the populations of the excited states due 
to cross-relaxation, the re-broadening will only occur until the 
radiative flux is large enough that the rate of stimulated emis- 
sion exceeds the rate for cross-relaxation. Additionally, in the 
case of the 22 GHz rotational H2O maser transition (616 - 523), 
the interaction between the different hyperfine components 
(F = 7 - 6, 6 - 5, 5 - 4, 5 - 6, 6 - 6 and 5-5) complicates 
a straightforward line analysis. Several papers have addressed 
this issue (e.g. Nedoluha & Watson 1991). The manifestation 
of the hyperfine components in the observed maser spectra also 
depends on the maser thermal line width, as well as the satu- 
ration and beaming angle. This is due to the fact that the hy- 
perfine components cause significant deviations from Gaussian 
symmetry. The deviations can be classified by a skewness pa- 
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rameter In this paper we examine the feasibility of using the 
skewness of the profile together with the line widths to disen- 
tangle the thermal line width of the 22 GHz masing medium 
and the level of saturation and beaming angle. We examine a 
sample of late-type stars and compare their high resolution to- 
tal intensity spectra to our non-LTE models. 

The H2O masers around late-type stars are thought to oc- 
cur in the expanding wind at a few hundred AU from the central 
star (e.g. Lane et al. 1987). The temperature of the masing gas 
is expected to be less than 1000 K, according to the analysis of 
H2O maser pumping schemes by Neufeld & Melnick (1990). 
As a result the intrinsic, thermal line width of the maser should 
be ^ 1 .5 km s"' . Larger line widths can be caused by turbulence 
of the masing gas. Model estimates of the intrinsic thermal line 
widths can thus give an upper limit to the temperature in the 
maser region and to the effect of turbulence. At a few hundred 
AU from the central star the outflowing material is still being 
accelerated. Therefore, as shown in Rosen et al. (1978), the ve- 
locity coherent paths in the tangential and in the radial direction 
are of similar length. As a result, for the H2O masers neither 
radially nor tangentially beaming necessarily dominates. The 
H2O maser spots have typical sizes of 10'^ cm (Reid & Moran 
1981), but as determined from VLBI observations they can be 
as small as a few times 10'^ cm (Imai et al. 1997). Due to beam- 
ing, the actual size of the H2O maser cloud is generally much 
larger than the observed size. Goldreich & Keeley (1972) have 
shown that for H2O masers a factor of 50 difference is likely. 
The maser path length could thus be as long as lO'"* - 10'^ cm. 
Along this path the velocity can change by a significant frac- 
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Fig. 1. Scaled intensity lines profiles for the non-LTE models 
with Vth - 0.6 km s"' (left) and Vth = 1.0 km s"' (right). From 
inside out, the dashed fine has emerging maser flux TbAQ. = 
10^ the short-dashed Une has TbAO = lO^" and the solid line 
hasTbAQ = lO". 

tion of the intrinsic thermal line width, especially near the edge 
of the envelope where the velocity gradient is steep. Here we 
also examine the efi'ect of such velocity changes on the shape 
of the 22 GHz H2O maser profiles. 

2. Background 

The rate equations for the populations of the upper and lower 
energy levels of the various hyperfine lines of the 22 GHz H2O 
maser have been solved along a uni-directional maser beam us- 
ing the method described in Nedoluha & Watson (1992, here- 
after NWa) and used in Vlemmings et al. (2002, hereafter V02). 
We have included only the three strongest hyperfine compo- 
nents (F = 7 - 6, 6 - 5 and 5 - 4), since it has previously been 
shown that the weakest of the hyperfine fines are negligible 
(Nedoluha & Watson 1991, Deguchi & Watson 1986). In V02, 
the rate equations are solved including a magnetic field, thus in- 
cluding all the different magnetic substates. The magnetic field 
does not noticeably influence the line shapes and line widths of 
the total intensity maser fines. Here we Umit ourselves to the 
case of no magnetic field. 

For the population levels of the upper states we solve 

= Mv) - (T + ry)niiv) + RijivXrijiv) - miv)) 

+m(^) f dvYgimiv). (1) 

A similar but reversed equation is solved for the population of 
the lower levels. 



Fig. 2. Line width as a function of emerging maser flux TbAQ 
for models with different thermal width. The thick solid line is 
the model for Vth - 1.0 km s~' with [F + Ty\= 5.0, the other 
lines are for [F + Fy]= 1.0. 

Here ni(v) and «/v) are the population numbers for respec- 
tively the upper (616) and the lower (523) rotational levels. / 
denotes the upper level hyperfine states {F -1 ,6 or 5), j is the 
related lower hyperfine state {F' - 6,5 or 4). The population 
levels are solved as a function of molecular velocity v and at 
different positions along the maser path. The statistical weights 
are designated by gi and gj, and are from Kukohch (1969). The 
pump rate /l,(v) is assumed to be the same for the different hy- 
perfine components and has a Maxwellian distribution. The re- 
sults depend on the the ratio A/1//1,- where Ai is the difference 
between the pump rate into the upper and lower levels. This ra- 
tio is of the order of a few percent (Anderson & Watson 1993). 
The rate for stimulated emission Rij is calculated using the lo- 
cal maser intensity and the hyperfine interaction coefficients as 
described in NWa. F is the decay rate for the molecular exci- 
tations. The cross-relaxation rate Fy describes the reabsorbtion 
of previously emitted infrared pump photons trapped in opti- 
cally thick transitions of the maser system, which generate a 
newly excited molecule at random, within a Maxwellian ve- 
locity distribution (0(v)). This has the effect of postponing the 
rebroadening of the maser line due to maser saturation. The 
rate for redistribution in velocity is not exactly the same for 
the redistribution among the different hyperfine components. 
However, both are expected to be roughly equal to the inverse 
lifetime of the state due to emission of infrared radiation. In 
NWa it was shown that the results scale by [flux/(r + F^)]. 
Most of the calculations are performed for [F h- Fy] = 1 .0 s"' al- 
though we also investigate higher values of F and Fy. A similar 
effect as the cross-relaxation due to reabsorbtion is produced 
by elastic colHsions between the maser molecules and inter- 
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Fig. 3. Skewness parameter as a function of emerging maser 
brightness TbAQ. The Hnes are labeled for different thermal 
width. The shaded areas indicate the 3cr errors when including 
observational errors as described in the text. 

mixed H2 molecules, that change the velocity but not the en- 
ergy of the colliding particles. This effect is described in detail 
in Elitzur (1990), where it was shown that due to elastic colli- 
sions, the narrowing of the maser line can be extended signifi- 
cantly. However, the analysis by Elitzur (1990) indicates that 
for maser transitions that have large infrared optical depths, 
such as the H2O maser discussed here, the effect of elastic col- 
lisions is small compared to that of F and Fy. 

Eq. [2 is iteratively solved along the maser beam, as the 
maser intensity itself, determined using the radiative trans- 
fer equations in NWa, depends on the level populations. We 
assume nearly one-dimensional maser propagation, with the 
beaming of the maser radiation represented by the solid an- 
gle AQ. The emerging maser fluxes will thus be represented 
in TbAQ, with Tb the brightness temperature. The beaming an- 
gle Afl is not a constant property throughout the entire maser 
regime. It depends on the maser optical depth |t| as beaming is 
more pronounced for the stronger masers. The beaming angle 
in the unsaturated cases is roughly proportional to |t| ', while 
in the saturated case it is proportional to |t|"^ (Elitzur 1992). 
All calculations were performed with TbAfl = 0.1 K sr, for the 
radiation incident onto the maser region. It was verified in NWa 
and V02 that the results are insensitive to the chosen value. 

2A. Profile shapes 

In Fig. [2 we show an example of the line profiles resulting 
from our calculations. The example profiles have been cal- 
culated for two dififerent values of the FWHM thermal width 
(vth) of the MaxweUian particle distribution. Assuming a ki- 



netic temperature T in Kelvin, v'th ~ 0.5 x (r/lOO)'''^ (NWa). 
Line profiles are presented for different values of emerging flux 
(TbAQ. = 10^10'° and 10"). Due to the hyperfine interac- 
tion the maser profiles show clear deviations from Gaussian 
symmetry with increasing maser brightness. The effects are 
largest for the smallest thermal widths. We also notice the re- 
broadening of the profiles when the maser starts to get satu- 
rated. This can also be seen in Fig.|2] which shows the maser 
line width as a function of emerging flux. In the unsaturated 
case the maser gets increasingly narrow, when the masers starts 
to get saturated the line width increases. The figure shows the 
relations for [F + Fy] = 1 and 5 s"'. As expected, the maser 
is able to remain unsaturated longer, for higher values of the 
decay and cross-relaxation rate. 

Aside from the narrowing and re-broadening of the maser 
lines. Fig. ^shows that due to the contributions of the hyper- 
fine components, the shape of the 22 GHz H2O maser profile 
is also a good indicator of the saturation level. The hyperfine 
components will cause the profile to be skewed toward the neg- 
ative velocities and the skewness will depend on the level of 
saturation. The deviation from Gaussian shape can be quanti- 
fied by the skewness parameter defined as: 

Skewness = - dv[ifr(v)- -]\ (2) 

^ J cr 

Here i^(v) is the profile of the maser line, cr the standard de- 
viation and vo the position of maximum intensity. The parame- 
ter is scaled by 'I', which is the area under the profile. A positive 
value of the skewness parameter signifies a distribution with an 
asymmetric tail out toward positive velocity, while a negative 
skewness parameter indicates a negative tail. Fig.^shows the 
skewness of the maser profile as a function of emerging bright- 
ness for different values of Vth, as determined from our mod- 
els. For the highly unsaturated masers the skewness parameter 
is strongly negative due to the multiple hyperfine components. 
The profiles become more symmetric when the ¥-1-6 hyper- 
fine component starts to dominate. When the maser starts to 
saturate, at TbAQ ^ 10"' for [F -1- F^] = 1 s"', the skewness 
parameter decreases again, because the ¥-1-6 hyperfine com- 
ponent is the first to saturate. As a result, the profile becomes 
more asymmetric. This is the most pronounced for the lowest 
values of v'th. For v'th < 1.0 km s"' the skewness parameter in- 
creases again for TbAQ > 10" after a strong decrease. For the 
higher thermal velocity widths the skewness parameter keeps 
decreasing slightly in the saturated regime. 

Because the skewness parameter is the third moment of the 
maser profile, it is strongly affected by observational errors. To 
include the influence of the errors, we have simulated our ob- 
servation by adding Gaussian distributed errors on top of the 
model profiles. We have used a signal-to-noise of 100, typi- 
cal for the weakest of the observed maser sources. The result 
is indicated with the shaded areas in Fig. |3l these correspond 
to 3cr errors. The effect of the errors is strongest in the un- 
saturated regime and decreases when the maser saturates. The 
spectral resolution of the maser profiles also affects the eiTors in 
the skewness determination. Higher spectral resolution obser- 
vations have smaller errors because of a more accurate deter- 
mination of the scaling parameter and velocity vo of the peak 
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Fig. 4. Skewness and line widths for from Vth = 0.8 (left) and 
1.0 km (right) when including a velocity shift AVm of 0.4 
(thin solid), 0.6 (short-dashed), 0.8 (dashed-dotted) and 1.2 
km s"' (long-dashed) as a function of emerging maser bright- 
ness TbAfi. The thick solid line indicates the result with con- 
stant velocity along the maser path. Due to the computational 
requirements we have not calculated the full brightness temper- 
ature range for the models which have a velocity gradient. 

intensity. Our models were determined with a spectral resolu- 
tion of 0.027 km s ' , similar to the observations. 

2.2. Velocity gradients 

Obviously the line profile can also be determined by velocity 
gradients along the radiation path, giving rise to observable 
skewness. Therefore, the calculations above have also been 
performed including a shift of velocity along the maser path 
(AVm)- This has been done by replacing the pump rate A{v) in 
Eq.n]for all hyperfine components by A{v + 6v), where 6v in- 
creases for every step along the maser path. The effects have 
been calculated for the models with Vth - 0.8 and 1.0 km s '. 
Results for the line widths and skewness parameter are shown 
in Fig.|4]for AV^ of up to 1.2 km s"'. In our calculations, the 
linear velocity gradient along the maser path is characterized 
by AVm, the velocity shift between the start and the end of the 
path. The small variations on the lines are due to the spectral 
resolution of our models, which is equal to the resolution of the 
observations. 

The skewness is very sensitive to AVm, especially in the 
unsaturated case. The effect of AVm becomes less for models 
with a relatively large thermal width. When AV^ becomes sig- 
nificantly larger than the thermal width, the skewness no longer 
depends on the intrinsic width. This explains that, in Fig.|3 the 
model for Vth - 0.8 and a AV^ = 1.2 km s ' is less positively 



skewed than the model with AVm - 0.8 km s"', while this is 
not seen in the models with Vth - 1 .0. Most models show less 
deviations of Gaussian symmetry in the saturated case. All the 
models in Fig. |3 were calculated with a negative velocity gra- 
dient along the maser path. A gradient in the positive velocity 
direction increases the negative skewness of the profiles for all 
cases. Since the front side of the envelope has the most nega- 
tive velocity, the velocity gradient through the envelope toward 
the observer is predominantly negative, for features originating 
from both sides of the circumstellar envelope. 

We find that a velocity gradient along the maser path in- 
creases the maser line width with an approximately equal 
amount in both the unsaturated as well as the saturated case. 
As soon as AVm approaches the intrinsic thermal line width, 
the measured line widths increase rapidly. Nedoluha & Watson 
(1988) have shown that for AVm a few times Vth, the maser line 
splits into several distinct, narrow features, each becoming an 
independent maser This has been observed for H2O masers in 
star-forming regions (e.g. Genzel et al. 1979, Walker 1984). 
Our observations of the circumstellar H2O masers do not show 
this line splitting, indicating that the largest velocity shift along 
the maser path is not more than ^2.0 km s"' . 

3. Observations 

The observations were performed at the NRAO' Very Long 
BaseHne Array (VLBA) at December 13th 1998 as part of the 
project designed to measure the circular polarization of circum- 
stellar H2O masers. The results of the polarization measure- 
ments are presented in V02. At the frequency of the 616 - 523 
rotational transition of H2O, 22.235 GHz, the average beam 
width is X! 0.5 x 0.5 mas. This allows us to resolve the dif- 
ferent H2O maser features in the CSE. As explained in V02, 
the data were correlated twice, once with modest (7.8 kHz 
= 0.1 km s"') spectral resolution, and once with high spec- 
tral resolution (1 .95 kHz = 0.027 km s"'). In this paper we use 
the total intensity profiles obtained with the high spectral res- 
olution data. We have performed 6 hours of observations per 
source-calibrator pair The calibrator was observed for 11/2 
hours in a number of scans equally distributed over the 6 hours. 
We used 2 filters of 1 MHz width, which were overlapped to 
get a velocity coverage of 22 km s"'. This covers most of the 
velocity range of the H2O maser 

3.1. Sample 

We observed 4 late type stars, the supergiants S Per, VY CMa 
and NML Cyg and the Mira variable star U Her. They are listed 
in Table. ^with type, position, distance and period. The sources 
were selected on the basis of 2 criteria; strong H2O masers and 
the availability of SiO and/or OH maser polarization observa- 
tions. The peak fluxes measured in the observations are shown 
in Table. ^ On the total intensity maps the noise is between 
~ 0.08 - 0.3 Jy. 

' The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agreement 
by Associated Universities, Inc. 
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Table 1. Star Sample 



Name 


Type 


RA 


Dec 


Distance 


Period 


V,ad 


Peak flux 








C ' ") 


(pc) 


(days) 


(km s-') 


(Jy) 


SPer 


Supergiant 


02 22 51.72 


+58 35 11.4 


1610 


822 


-38.1 


76.1 


UHer 


Mira 


16 25 47.4713 


+ 18 53 32.867 


277 


406 


-14.5 


12.5 


NML Cyg 


Supergiant 


20 46 25.7 


+40 06 56 


1220 


940 


0.0 


48.2 


VY CMa 


Supergiant 


07 22 58.3315 


-25 46 03.174 


1500 


2000 


22.0 


244.1 



Table 2, Results 
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15.0 


17.8 


0.94 


0.01 


+0.02 


-0.05 




18.3 


58.5 


0.80 


-0.36 


+0.06 


-0.08 




24.6 


126.2 


0.60 


0.07 


+0.02 


-0.03 




25.2 


65.6 


0.79 


0.34 


+0.07 


-0.06 




27.1 


25.3 


0.88 


0.22 


+0.06 


-0.06 




28.8 


43.1 


0.61 


-0.32 


+0.05 


-0.06 


NML Cyg 


-20.0 


3.66 


0.54 


0.64 


+0.32 


-0.36 




-20.2 


5.42 


0.54 


-0.14 


+0.11 


-0.30 




-20.4 


10.5 


0.40 


-0.09 


+0.08 


-0.19 




-21.2 


48.2 


0.66 


-0.03 


+0.02 


-0.03 


UHer 


-17.6 


2.91 


0.50 


5.47 


+ 1.65 


-1.60 




-17.6 


1.16 


0.73 


0.69 


+0.55 


-0.63 




-17.8 


12.20 


0.75 


0.69 


+0.18 


-0.16 




-19.2 


2.46 


0.47 


-0.63 


+0.38 


-0.93 




-19.2 


2.07 


0.51 


-0.14 


+0.33 


-0.76 




-19.3 


1.36 


0.51 


-0.27 


+0.46 


-1.14 



* Bad spectral channels affect skewness determination 



4. Results 

The results of the observations are shown in Table|3 The table 
gives the LSR velocity of the maser feature, the maser strength 
in Jy and the line width and skewness parameter. The features 
labeled * have a few bad spectral channels in the wings of the 
profile due to interference. This can have a significant effect 
on the skewness of the profile; we have verified that the val- 
ues listed are the upper limits for these features and the actual 
skewness could be up to a; 1 .0 lower. The quoted 3cr error val- 
ues of the skewness, are based on the errors in the determina- 
tion of v'o and cr for use in Eq.|2| These depend on the signal-to- 



noise ratio of the maser feature, and on the spectral resolution 
of the observations. 

We will plot the observed line width and skewness as a 
function of brightness temperature in order to make a com- 
parison with our models. One should realize that there is not 
necessarily a single value for width and velocity shift that is 
valid for all maser features in any star. It is natural to expect 
the results to be quite scattered. 

Figs. |5l and |6l show the observed line widths Av^ and 
Skewness parameters plotted against the logarithm of the maser 
brightness temperature T\,. This is calculated using the fluxes 
hsted in Table|2using the beam size of the observations as the 
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Fig. 5. Line widths of the observed maser features as a function 
of emerging flux. The error bars are the 3cr errors. 
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Fig. 6. Skewness parameter of the observed maser features as 
a function of emerging flux. The error bars indicate 3cr errors 
and the dashed-dotted line is the weighted average. 



maser spots are unresolved, thus T\, is formally only a lower 
limit. However, at the distance of the star the beam has a diam- 
eter of approximately 1-2 x lO'"^ cm and is comparable to the 
typical maser spot size (Reid & Moran 1981). We thus expect 



the determined brightness temperatures to be with the real 
values. 



4.1. UHer 

For the Mira variable star U Her most of the maser features 
show line widths of w 0.5 km s ' with 2 of the features show- 
ing widths of X 0.75 km s '. There is no clear trend with 
increasing maser brightness. The H2O masers around U Her 
have also been observed with the Very Large AiTay (VLA) by 
Colomer et al. (2000). They observe line widths several times 
wider than found in our observations, but the VLA is unable to 
resolve the individual spots and they only used a spectral reso- 
lution of 0.33 km s ' . It is likely that the larger line widths are 
due to the blending of several smaller components. 

Because U Her has the weakest masers, the eiTors on the 
skewness parameter are relatively large. The weighted aver- 
age of the skewness is 0.55, a possible indication of large 
AVm. The skewness is positive for the maser features with 
Vlsr ~ -17.5 km s"', while it is negative for the features with 
Vlsr ~ -19.0 km s ^ All the maser features are blue-shifted 
with respect to the stellar velocity of Vlsr = -14.4 km s'^. 
As discussed above, AVm will be negative because the maser 
exists in a radially accelerating outflow. This results in a fre- 
quency shift along the maser feature in the direction of the hy- 
perfine components. If the features with a velocity closer to the 
stellar velocity are located closer to the star, they will experi- 
ence more acceleration along the maser path and thus a larger 
velocity gradient. This can explain that the skewness of these 
features is more positive than the that of the other features. The 
maser feature with the skewness of 5.47 however, is likely a 
blend of two overlapping maser features. 

4.2. NML Cyg 

The line widths of the maser features around NML Cyg are 
found to lie between 0.4 and 0.7 km s"^ Again no clear trend 
can be determined, especially since we have only been able to 
detect 4 maser spots. MERLIN observations by Richards et al. 
(1996) show a much more complex structure with many fea- 
tures. The line widths they find cover a large range between 0.4 
and 2.5 km s"'. 

The weighted average of the skewness is -0.01. The line 
widths and skewness, indicate an intrinsic thermal width of 
between 0.8 and 1.0 km s"' and AV^ of ^ 0.75 km s"^ A 
larger AVm would result in a more positive skewness, while a 
larger thermal width does not allow for the observed narrow 
line widths. 



4.3. S Per 

The line widths of the maser features around S Per show 
some indication of an increase with maser brightness. The line 
width increases from 0.4 to 0.9 and varies approximately as 
Av oc {Tbf'-^. However, 2 features show a width which is higher 
than expected. But, as seen in Fig. |3 this can be due to a 
velocity shift along the path which is more than 1.2 km s 
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Alternatively, the actual spot size of these features could be 
less than 10'^ cm leading to an underestimation of T\,. The 
skewness of the other features indicate a AV^ of not more 
than 0.8 km s"'. Similar velocity shifts have been observed by 
Richards et al. (1999). The observed line widths correspond to 
a thermal line width between 0.8 and 1 .0 km s ' . 

A combination of both the skewness and line width analy- 
sis allows us to estimate, from our models, the actual emerging 
maser brightness temperature and beaming angle. For the es- 
timated velocity shifts along the maser path and thermal line 
widths, our models give observed line widths of 0.4 - 0.9 
km s"^ for T^AQ. between 10'° and 10" K. With observed 
brightness temperature T\, - lO'^-lO'^, this indicates that the 
beaming solid angle AQ ^ 10"^ sr. 

4.4. VYCMa 

The line widths of the H2O masers around VY CMa show a 
clear dependence on the brightness temperature. The line width 
varies with Av oc (Tb)"", where a ~ 0.3. A similar decrease 
of line width with intensity was observed before on the maser 
flare in Cepheus A (Rowland & Cohen 1986), though there a « 
0.5. We observe Hne widths from 0.6 to 0.94 km s ' . MERLIN 
observations (Richards et al. 1999) have found many features to 
have a velocity extent of up to 3 km s ' , indicating the presence 
of large velocity shifts along the maser. 

For VY CMa, the skewness and line widths indicate a ther- 
mal width in the H2O maser region of x 1.0-1.2 km s 
with AVm up to 1 .25 km s" ' . Using this, according to our mod- 
els, T\,AQ. is between 10^° to 10'^ K, while the observed lower 
limit for ^ 10'^ K. This indicates a beaming solid angle of 
* 10-^ sr. 



5. Discussion 

The non-LTE method used here to model the shape of the 
22 GHz H2O maser spectrum was previously used to deter- 
mine the circular polarization of the maser in the presence of 
a magnetic field (V02). There the non-LTE modeling method, 
is compared with a standard LTE method which is described 
in Vlemmings et al. (2001). The non-LTE approach was devel- 
oped to analyze circular polarization spectra and successfully 
reproduced the distribution over magnetic components. In prin- 
ciple, similar physical processes are important for understand- 
ing the line shapes of the H2O maser blend and we apply the 
same analysis here. Recently, Watson et al. (2002) have used a 
similar method to examine the line shapes of interstellar H2O 
masers, although their analysis has been based on only a two- 
level maser transition, thus only 1 hyperfine component is in- 
cluded. Instead of using the skewness of the maser profile, they 
examine the second order deviations of Gaussian symmetry, 
the 'kurtosis'. They find that the interstellar H2O masers are 
found in gas with temperatures greater than 1200 K, for which 
the influence of multiple hyperfine components is thought to 
be small. Since the gas in CSEs is at a lower temperature (e.g. 
Goldreich & Scoville 1976), the efi'ect of the other hyperfine 
components is larger, as is indeed seen our results. 



The large spread in the observed skewness parameter indi- 
cates that velocity gradients along the maser path play a large 
role in determining the shape of the H2O maser line profiles. 
35% of the maser features have a positive skewness, which only 
occurs for masers that are starting to saturate and which have a 
relatively large negative AVm. 

The observed Av oc {Ti,y dependence, with different val- 
ues ofa(a^ -0.3 for S Per and 0.3 for VY CMa), can be ex- 
plained by the relation between beaming angle AQ. and maser 
optical depth |t|. As noted before, in the completely unsaturated 
case AQ oc |t|"', while this increases to AQ oc |t|"^ in the satu- 
rated case. Additionally, in the unsaturated regime the amplifi- 
cation is exponential, so T\, oc e^'"!, while in the saturated regime 
the amplification becomes linear and oc |t|. This means 
that, while unsaturated, the beaming angle AQ oc (In Tb)"', 
and when saturated AQ oc T^^. Meanwhile, our models in 
Fig. 12 show, that before the maser starts re-broadening the 
line width is almost constant with increasing TbAQ. after re- 
broadening, when saturation occurs, the line widths becomes 
approximately proportional to (TbAQ)"-^. Using the results for 
the variation of AQ with Th we then find that for unsaturated 
masers Av oc (Tb)"'^, while the line width of saturated masers 
will be proportional to (Tb) The value of a can thus in- 
crease from X -0.7 in the mostly unsaturated regime to 0.5 
in the saturated regime. This seems to indicate that the masers 
around S Per are still mostly unsaturated, while those around 
VY CMa are close to becoming saturated. However, the pres- 
ence of velocity shifts along the maser path complicates this 
straightforward explanation. For example, if the high bright- 
ness maser features trace the lower velocity shifts because a 
strong gradient disrupts the maser growth, we can also ob- 
serve a decrease of the line width with increasing brightness. 
Additionally, velocity shifts will introduce a big spread in the 
relation between line width and maser brightness. 

An additional uncertainty in the analysis is introduced by 
assuming a linear maser model. While astrophysical masers, 
and especially H2O masers, are thought to be highly elongated, 
the actual three-dimensional shape of the maser will affect the 
maser line width and shape as described by Elitzur (1994). 
There it was shown that the linear maser model is sufficient 
to describe the maser profile out to ~ 2vth for spherical and 
~ 3-5vth for elongated cylindrical masers. However, crossing 
maser beams will suppress the emission in the outer wings of 
the maser profile. This can cause naiTowing of the maser line 
to be observed even for saturated masers. In the case of the ob- 
jects studied in this paper, the effect is expected to be small or 
even negligible, as the suppression of the maser profile wings 
only becomes significant when the maser becomes saturated. 
Then, as shown in Elitzur (1994), the suppression effects on the 
profile are only noticeable when the dynamic range of the spec- 
trum is large. For a filamentary maser, the dynamic range has 
to be higher than ^ 1 x 10^, while the largest dynamic range for 
our maser sources is w 7 x 10^ for the strongest maser feature 
of VY CMa. If the masers were purely spherical however, the 
measured linewidth decrease for increasing maser brightness 
for the VY CMa H2O masers could be partly due to this geo- 
metric effect. As a result the estimated intrinsic thermal width 
and the beaming solid angle might be slightly underestimated. 
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The presence of large velocity gradients along the maser 
path, and the uncertainty in the exact maser geometry, increases 
the difficulty of determining underlying thermal widths, level 
of saturation and beaming angle. However, from comparison 
with the models we have been able to give estimates for most 
of the sources presented here. In general the thermal line width 
of the maser medium in our sources is between 0.8 and 1.0 
km s~^ . This line width can be caused by the actual temperature 
of the maser gas, as well as by turbulent motions. The H2O 
maser region should not have temperatures higher than T ^ 
1000 K (Neufeld & Mehiick 1990), which impUes that the Une 
widths due to the thermal motions should not be higher than 
» 1.5 km s"^. This is consistent with our observations, which 
indicate a temperature of 7 « 400 - 600 K. The Une width 
due to turbulent motions should then be of the same order of 
magnitude or less. This is significantly less than the turbulent 
width determined by Colomer et al.(2000), who give a width of 
the order of 2 - 4 km s~'. But as discussed above they found 
larger line widths, possibly due to the fact that the VLA was 
not able to resolve individual maser features. 

Similar to earlier observations (e.g. Spencer et al. 1979), 
our observations suggest that the circumstellar H2O masers 
are either unsaturated or only slightly saturated, with intrin- 
sic brightness temperatures between - 10'" - lO" K. From 
the models we can estimate the typical beaming angle to be be- 
tween AQ = 10"'^ - 10"^ sr, similar to values found before (e.g. 
Nedoluha& Watson 1991). 

6. Conclusions 

A detailed analysis of the line profiles of 22 GHz H2O masers 
can be a valuable tool to estimate saturation, thermal veloc- 
ity width and velocity shift along the maser path in the maser 
region. Because the 22 GHz transition consists of multiple hy- 
perfine components, the maser profile shows clear deviations 
from symmetry at various stages of saturation. Unfortunately, 
the efi'ect of velocity shifts along the maser path, can not eas- 
ily be detached from the effects of hyperfine contributions. 
Additionally, for VY CMa, the actual unknown maser geom- 
etry introduces another level of complexity. 

However, we have been able to obtain some estimates on 
the thermal line widths and beaming angles on the H2O masers 
around a sample of 4 late-type stars. It is clear that velocity 
shifts along the maser path of up to 1.25 km s~' play an im- 
portant role, and will have to be taken into account for further 
comparisons between models and observations. The thermal 
line widths in the maser regions are found to be mostly around 
1.0 km s"\ indicating temperatures of approximately 400 K in 
the masing region. Most masers are thought to be only slightly 
saturated at the most, with TbAQ % lO'^ K. Beaming solid an- 
gles for the circumstellar H2O masers are found to be typically 
between 10"-^ and 10"^ sr. 
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